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Abstract 

As a nuclear spin model of scalable quantum register, the one- dimensional chain of the magnetic atoms 
with nuclear spins 1 /2 substituting the basic atoms in the plate of nuclear spin free easy-axis 3D antiferromag- 
net is considered. It is formulated the generalized antiferromagnet Hamiltonian in spin-wave approximation 
(low temperatures) considering the inhomogeneous external magnetic field, which is directed along the easy 
axis normally to plane of the plate and has a constant gradient along the nuclear spin chain. Assuming a 
weak gradient, the asymptotic expression for coefficients of unitary transformations to the diagonal form 
of antiferromagnet Hamiltonian is found. With this result the expression for indirect interspin coupling, 
which is due to hyperfine nuclear electron coupling in atoms and the virtual spin wave propagation in 
antiferromagnet ground state, was evaluated. 

It is shown that the inhomogeneous magnetic field essentially modifies the characteristics of indirect 
interspin coupling. The indirect interaction essentially grows and even oscillates in relation to the interspin 
distance when the local field value in the middle point of two considered nuclear spin is close to the critical 
field for quantum phase transition of spin-fiop type in bulk antiferromagnet or close to antiferromagnetic 
resonance. Thus, the external magnetic field, its gradient, microwave frequency and power can play the role 
of control parameters for qubit states. 

Finally, the one and two qubit states decoherence and longitudinal relaxation rate are caused by the 
interaction of nuclear spins with virtual spin waves in antiferromagnet ground state are calculated. 

Keywords: Decoherence, easy-axis antiferromagnet, indirect coupling, inhomogeneous magnetic field, nu- 
clear spin, quantum register and qubit. 
PACS numbers: 75.10.Pq, 82.56.-b, 75.60.-k, 75.50.Ee 

1 Introduction 

In 1998 B.Kane proposed the scheme of large-scale (scalable) nuclear magnetic resonance (NMR) quantum 
computer with the register in the form of a regular chain of atoms "^^P with nuclear spins 1/2 as qubits 
implanted into a near-surface layer of nuclear spin free substrate ^*Si (Ref.|;lJ). It was assumed that the control 
of individual nuclear spin-qubits states is performed by electrical gates A, which change the electron- nucleus 
hyperfine interaction in atoms '^^P. The indirect nuclear spin coupling for neighboring donors is controlled by 
gates J, which change the overlap degree of donors wave functions (exchange integral). Two-qubit operations for 
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Figure 1: The Kane's scheme of quantum register on nuclear spins of donor atoms ^^P. S is electron spin, I is 
nuclear spin (shot arrows), B is external magnetic field, V is electrical potential on gates J. The gates A are 
not shown. 
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Figure 2: The scheme of nuclear spin quantum register based on one dimension antiferromagnetic chain of 
two-lag ladder type [6] 

far spatially separated qubits can be produced using SWAP quantum operations between neighboring qubits. 
The separation between neighboring donors in this scheme should be no more than 20 nm (Fig.[T]). 

The main difficulties encountered in realization of Kane's scheme were connected mainly with the necessity 
of using very low spin temperatures for nuclear spins, using high precise technology for realization of quantum 
register structure and the complexity of quantum operation performing on states with random access to any 
qubits. As an alternative to this semiconductor scheme, a model for antiferromagnet-based NMR quantum 
cellular automata, which excludes the same difficulties of Kane's scheme, was proposed in our papers [2I3|. 

Authors of paper [Ij have proposed a variant of scalable quantum computer model, which does not use 
the gates of type J. For quantum-information exchange between electron states of removed atoms "^^P the 
implementation of spin wave propagation in the thin antiferromagnetic layer with easy axis of anisotropy grown 
on silicon layer with implanted donor atoms '^^P, was suggested. For the local generation of spin wave at input 
of register and for signal detection in readout procedure the use of the scanning tunneling microscope was 
suggested. 

Another model of NMR quantum register based on two-leg ladder ID antiferromagnet chain, where nuclear 
spins-qubits are placed in the inhomogeneous magnetic field and separated by a distance up to several tens of 
lattice constants, was proposed in papers [Sl[51[7] (Fig. [2). Relatively large inter-qubit distance in this case is 
effective both to diminish the nuclear dipole coupling between qubits and to give the considerable decoherence 
suppression. It is assumed, that to achieve the required difference of resonance frequencies needed for individual 
access to individual qubits the gradient of external magnetic field along the spin chain can be used. It is suggested 
also that two qubit operations on this register may be performed by switching the indirect interaction between 
longitudinal components of removed nuclear spins-qubits. To do this, it is assumed to use of excitation of the 
spin wave packets by microwave pulses. The authors consider as the possible candidates for antiferromagnetic 
base a number of organic materials having energy gap for spin excitations. They should consist of atoms of 
stable nuclear spin free isotopes hke (98.90%), i^-i^O (99. 962%), (99.634%) (in brackets their natural 
abundance is given). This also may be in particular same Haldane systems (Ref.[6l [7]). 

At very low temperatures and in the absence of external magnon excitations the indirect interaction between 
longitudinal components of nuclear spins does not occur. However, the transverse part of interaction is different 
from zero and is determined by hyperfine nuclear spin and virtual magnon interaction. The theory of indirect 
nuclear interspin coupling in homogeneous isotropic antiferromagnet was the first given by T.Nakamura (Ref.[i8i]) 
and than developed also in Refs.[9l fTOl ITT]. 

We have considered a model of NMR quantum register in the form of one-dimension chain of atoms with 
nuclear spins / = 1/2, placed on the thin plate of nuclear spin free 3D antiferromagnet with easy-axis anisotropy 
at low temperature. It is proposed that the external field is directed along the easy axis, which is normal to the 
plane of the antiferromagnet plate. In the case of homogeneous external magnetic field, the energy spectrum 
of spin waves (magnons) in easy axes antiferromagnet for field B lower than the same critical value Be is 
characterized by the energy gap. We have shown (Ref. T7) that close to the critical point of spin-flop quantum 
phase transition the range of indirect nuclear spin coupling can run up to a great value. In the next papers 
[Tni[TH[TS], the previous model was extended into the case of external magnetic fleld with gradient along nuclear 
spin chain. 

In present paper, we have presented the further investigations and development of this model, in which, 
as distinct from paper [13j, the role of umklapp processes between distinct cells of reciprocal sublattice were 
accounted. Corrected and more refined analytical expressions were obtained and numerical calculations for 
the interspin indirect coupling were performed. It was shown that the character and magnitude of indirect 
coupling between nuclear spins essentially changes as the local external field in its middle point comes close to 
the critical field for homogeneous phase transition. The indirect interspin interaction at the antiferromagnetic 
resonance condition was also calculated. As a result, we have new fresh possibilities for control not only the 
individual nuclear spin resonance frequency, but also the interaction between separated spins without resort to 
many gate-controlling systems. The one-qubit and two-qubit decoherence and longitudinal relaxation rates are 
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Figure 3: The scheme of antiferromagnet based nuclear spin quantum register in external field lower than for 
spin-flop phase transition B (x) <C Be- The oriented (that is nonprecessing) arrows represent here the ground 
states of corresponding individual Bloch vectors. The nuclear spins Ia (shot red arrows) are contained here 
only in atoms of sublattice A 

caused by the interaction of nuclear spins with virtual spin waves in antiferromagnet ground state were finally 
considered. 

The paper is organized as follows. In Section 2 we discuss the structure of considered model with the 
Hamiltonian of easy axis antiferromagnet plate in inhomogeneous magnetic field. In Section 3 we formulate 
the generalized antiferromagnet Hamiltonian in spin-wave approximation (low temperatures) with a constant 
gradient along the nuclear spin chain. In Section 4 the diagonalization of antiferromagnet Hamiltonian is 
carried out. In Section 5 asymptotic expressions for coefficients of unitary transformations to the diagonal form 
of antiferromagnet Hamiltonian are found for small gradient parameter. In Section 6 we evaluate the expression 
for indirect coupling between transverse components of nuclear spins. In Section 7 the indirect interaction, close 
to antiferromagnetic resonance condition, is evaluated. In Sections 8 and 9 the calculation of one and two qubit 
states decoherence and longitudinal relaxation rate are presented. 

In Conclusion 10 we discuss the some prospects of considered quantum register model. Some details of 
calculation are placed in the Appendixes. 

2 Spin Hamiltonian of easy-axis antiferromagnet plate in inhomo- 
geneous magnetic field 

The antiferromagnet electron spin system of considered model plays here the role of an environment for nuclear 
spin quantum register, whose interaction with spin wave leads on the one hand, to indirect coupling between 
nuclear spins and on the other hand to decoherence and relaxation processes of their states. 

It may be proposed to use the natural antiferromagnetic crystals with easy-axis anisotropy as an antiferro- 
magnet plate (or film, d is its thickness). As examples, they may be crystals CeC2 with tetragonal and FeCOa 
(siderite) with trigonal symmetry. The basic isotopes of these crystals ^^C, ^^Fe (91.7%), ^^O, "^'^^Ce (99.6%) 
have no nuclear spins (in brackets percent isotopic abundance is given). To form the one-dimension nuclear spin 
chains the isotopic substitution atoms, such as ^^C in corresponding crystal lattice sites, for isotopes ^'^C with 
nuclear spins 1/2, are proposed. One would expect that period of such solid state NMR quantum registers may 
be much more than periods of crystal lattice. 

The simple antiferromagnet model to be studied here consists of two incorporated to each other tetragonal 
magnetic sublattices A and B with N — N±Nz sites in each sublattice, where N± = N^Ny ^ 1 are the sites 
numbers in plane of plate (x,y-axea) and > 1 is the sites numbers in z direction. The atom sites of sublattices 
are numbered respectively by numbers j and i. Each sublattice constant in the plane of the plate is a± and 
along symmetry axis is Uz- 

The external magnetic field B (x) in considered model (Fig. ^ is directed parallel to z-axis and to the 
easy axis. As calculations (Ref.|16)) show, the field gradient dB^ {x) /dx = G of the order of 1.4T//im, may 
be obtained by using dysprosium micro magnet with dimensions 10 x 4 x 400 m'^ at a distance of 2.07 m. We 
assume here that the field has somewhat lesser constant gradient G ~ 0.1 T/ fim along nuclear spin chain (x-axis 
in plane of the plate) . This value of the gradient corresponds to the difference of resonance frequencies of the 
order of 100 kHz for two nuclear spins, being separated by 100 a_L {a± ~ Inm). 

The qubits number in quantum register will be limited by planar structure dimensions. For example, the 
structure with linear dimension of the order of 10 jim will have 100 qubit register with period L = lOOa^. 

The starting electron spin Hamiltonian of 3D easy-axis antiferromagnet with interaction only between neigh- 
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bouring atoms, which belong to the distinct sublattice, is represented for our model as 

(N ^ \ 

J2 B {x^) Saz (r,) + ^ S (x,) Sbz (r,) j + (1) 

N Z 

+27s?i/Z ^ ^ { Sa (r.) Sb (r, + r^) + Ba Saz (r,) ^b. (r^ + r^)}, 

i 5 

where Sa (fi) and Sb + ys) are electron spin operators {S = 1/2) for neighbouring sites of sublattices A and 
B, B {xi) = B + Gxi, B is the field value at the origin of the coordinates Xi, Z ^ 6 is the number of neighbouring 
atoms for tetragonal sublattice, js — 175.88 rad GHz/T (75/27r = 28GHz/T) is electron spin gyromagnetic 
ratio, h = 1.054 • lO'^^ Js/rad. 

Note, that the product of spin operators in Eq.([T]) written in matrix representation is common designated 
by symbol 0. In the following this symbol will be for brevity omitted. 

In Eq.© the parameters ~ 10 - 100 T, Ba ~ 10"^ - 1 T > 0, Be = ^/2BeBa + B\ are exchange 
field, anisotropy field and spin-flop field for easy-axis antiferromagnet (particularly, for FeCOs: Be = 35 T, 
Ba = 3.3 T, Be w 15.5 T). 



3 Spin Hamiltonian of antiferromagnetic plate in spin wave approx- 
imation 

The radius vectors of positions for magnetic atoms (for example, for atoms A) will be presented as — > 
{a±Xi, a±yi^,azZi), where {xi,yi),Zi ^ r± , z are now dimensionless radius-vector values. The dimensionless 
volume of elementary sublattice cell and of reciprocal sublattice cell in this case becomes equal to unity {va = 1) 
and the plate dimensions in three perpendicular directions are coincident correspondingly with the site numbers 

Because of large values N^, Ny, antiferromagnet spin states are essentially independent of boundary condi- 
tions at interface plate; it is convenient to use the Born-Karman periodic conditions for them. 

The local spin operators of distinct sublattice atoms are expressed in terms of wave representation as follows: 

Sax (r») - iSAy (rj) = 5^ (r^) = a+ = l/\/iV^ a+ exp ( iqr.i), 

q 

Sax (r*) + iSAy (rO = (r-O = ai ^ 1///V^aq exp (-^qr^), (2) 

q 

(r,) 5^ b, - 1/%/^^ 6q exp (-iqr,), S+ (r,) ^ b+ = 1/%/iV^ 6+ exp ( zqr,), 
q q 

with 

SAz{r,)= 1/2 -a+a,= 1/2 - (l/A^) ^ a+ exp (z (q' - q) r,) , (3) 

q,q' 

Sbz (r,) = -1/2 + b+b, = -1/2 + (1/iV) ^ 6+ 6q exp {i (q' - q) r,) . 

q,q' 

The operators Oq and Oq act as operators of creation and annihilation of spin waves (magnons) with di- 
mensionless wave vector q (magnon quasi-moment) for sublattice A and operators 6q and 6q accordingly for 
sublattice B. These operators do not satisfy the Bose commutation relations. However, it is known that by 
means of Holstein-Primakoff transformation they be converted to magnon operators with boselike commutation 
relations. 

The operators in site representation afa-i, and b'^'bj in the second terms in Eq.Q have the eigenvalues and 
1. They describe the number of flopped spins for one site of sublattice. The mean electron spin state deviations 
at sufficiently low temperature in vicinity of ground state T <^ 75^1 {Be — B) /k^ (k^ is Boltzmann constant) 
are smafl, that is (0|a+a,|0) < 1, (0|6+5i|0) < 1. 

In this case, one may belief that spin waves are created predominantly by the transverse components of 
electronic spins (therein lies the spin- wave approximation) . Spin wave operators take approximately the boselike 
form without using the Holstein-Primakoff transformation. The commutation relations for local Pauli operators 
take form (Ref.|I7j): 
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[ai, a'j;] ^ S^^i', [bi.b:^] = Si^i', (4) 

[ui, a.^,] = [a,;, bi>] = [ai,5+] = 0, [fei, a^] =0. 

Suggesting than that spin states modes have equal values at the interface of plate, we will have for transverse 
component of wave vector for thin plate the discrete values {d is the thickness of plate): 

(n) = 27rn,/iV„ 71, = 0, 1, 2, . . . , TV, - 1, TV, = d/a, > 1, (5) 

that is these spin modes in z-direction are the standing wave. Let us consider the plates with thickness d 
whereby the nonzero transverse modes with 7^ are not excited (conditions see Sec. 7) and the zero-order 
transverse spin mode with = is only populated. In this case, the spin waves are propagated only along the 
plate surface with two-dimensional wave vectors qx = {Qx^Qy), where q^-, qy = 2iTnx,y/Na;^y. For even N^^y we 
have values Ux^y — 0, ±1, ±2, . . . , ± {Nx,y/2 — 1) , Nx^y/2, that is components of wave vector sweep all states 
in one cell of reciprocal sublatticc (the first Brillouin zone) with < \qx\, \qy\ < tt. The sites of atoms in the 
sublattice A (and similarly in sublatticc B) are determined now by two-dimensional vectors rx,i = {xi, yi). 

Let us go from three-dimensional magnon operators to two-dimensional operators a+ = (^qj^^g.fij = 
b^^dq_,fi. As a result, we come to a simplified quasi- two-dimensional Hamiltonian for the antiferromagnetic 
plate. 

Then taking into account that the distribution of plane wave vectors for N^.y 3> 1 is practically contin- 
uous, we will go from the sums to integrals and from Kronecker's symbol to Dirac's (5— fimction (Appendix 2 in 
Ref.[I7]): 



q-L,92=o 



^(2^)ViV±'5(q^-q'^). 



In the general case, magnon quasi- moment q_L is defined modulo a vector of the reciprocal plane sublattice Q^. 
This essential property was ignored in our early works [121 113] . Considering only = and by summing over 
sites of direct sublattice in terms of Hamiltonian ([T]) after substitution of Eqs.([2]), ([3]), we will have: 

N N± 

1/N exp {i (q' - q) r.) = ^ exp {i (q'^ - q^) r^,,) = ^ '5(q'^_q,), , • (7) 

1=1 i=l Q 

where Qj. = {Qx = 2T:mx, Qy = 2TTmy), rux^y = 0, ± 1, ±2, . . . , ±cxd. 

Because the external magnetic field is inhomogeneous along x-axis, quasi- moment components qx^ unlike 
components g^, g,, are not integrals of motion and magnon states should be defined as superposition of 
states for qx and q'x — ~'r Qx- At the same time, displacement of components qy on Qy represents sim- 
ply umklapp processes of states in other elementary cells of reciprocal sublattice, which represents the iden- 
tical transformation and therefore then it will be assumed Qy — 0. Thus in the following, only the values 
Q-L = {Qx = 2nmx, Qy — 0) will be considered. Taking in to account the above-mentioned remarks we will 
write Eq.Q in the form 

E'^(ql-q^),Qx E '^(-3i-fc),QJ9;, 9«- (8) 

Let us go also in the Eq. ([8]) to continuous variables for qj^ but retain discrete values for Qx ■ We will obtain: 
]^,MQJi',,iy ^5](2^)VA^±'5(9i-fe-Q.)<5(9;-<Z,). (9) 



E' 



and 



X, exp {i (ql - qx) r^,) ^ -i (27r)' /N^ ^ dS {q'x - qx - Qx) /dqx 5 {q'y - qy) . (10) 



i=l 



We will then go to magnon operators with continuous arguments q^. For example, for creation operators 
of magnon state with wave vector q^ in the range from q^ to q^ -|- dq± , we will write 

^ 27r/^a+ (q^) , ^ 27r/^b+ (q^) . (11) 
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By the use of Eqs.©, © and (fTU|) . we will find the expressions of sums, in spin Hamiltonian ([T]) involved 
(in a like manner for operators 6^, , bqj^ and other operator pairs): 

i=lqj_,q'^ q±,Qx 

l^a^ {qx + Qx,qy) a (q±) dq±, 
-f -i jYa+{qx + Qx,qy)da{q±)/dqxdq±, (12) 



JVi 



(V^-l)E E <^^'tq,exp(z(ql+q^)r^,,)l/6 

i=l qj_,q'^ 

^ lYlq±a,^ iqx + Qx,qy)b+ {~-q±)dq± 
where for simple tetragonal lattice {Z — 6, rs — (±1/2, ±1/2, ±1/2)) 



Eexp {i {q±r±s)) + 2 
L s 



7qi =7q^ =7qi,g.=0 = 1/6 



Eexp {i {q±r±s)) + 2 

L <5 



(cos (g:r/2) + cos {qy/2) + 1) /3. 



After passing in Eq.([T]) to dimensionless designations: 

jsBe =loe, < Ba/Be = foA < 1, < B/Be = b, Gg^/Be = g^ IQ-^ , 



(13) 



(14) 



we will obtain the expression for spin Hamiltonian of considered easy-axis antiferromagnet model in the vicinity 
of ground state up to quadratic members relative to operators a (qj.) and b {q±) (spin- wave approximation) in 
the form: 

hs - Hs/ {Hlue) ^-N{1 + &a) /2+ 
+ jY.{{l + bA-b) a+ (q^ + Qx) a (q^) + (1 + + 6) b+ (-q^) b (- (q^ + Q^)) + 



^5 [a+ (qx + Qx) da (qx) /9fe - (96+ (-qx) /dqx) b {- (qx + Qx))] + 
+7q^ [a+(qx+Q±) 6+(-qx)+a(q±) 6 (- (qx + Q±))] } rfq±. 



(15) 



Hamiltonian (flSl) takes into account here the correlations of niagnon states for distinct cells of reciprocal 
sublattice, which are expressed as summation over all vector components Qx- In early works [T51 [TH] this 
circumstance was not considered. 

With the help of Eq. ^ , ([6]) and ([7| the commutation relations for operators a+, , a^^ take the form (similar 

to operators 6+^, , 6-q^): 



TV, 



Eexp(i(q'xrx,j - q^rx,,')) = E 



9i. <}x+Qx'^<Ji, 9j, J 



or 



E'^9i + Qi>«-+Q-'^9i^9!'- 

Going in Eq. (|17p to the continuous arguments qx, qy for magnon operators we will obtain: 
[a(qx),a+(gL + QL,9;)] = E S {q', - qx - Q x + Q'x) S {q'y ~ qy) , 



(16) 
(17) 

(18) 
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Q.-Q' 

Using the commutation relations (|18p and notation — = Q", we will obtain for the magnon operators 
in Heisenberg representation a (r, q) and (r, q) the following system of differential equations (here r = LUEt 
is dimensionless time): 

ida (r, q±) /dr = [a (r, q_L) , hs] ^ 
= Y,{i^ + bA-b + ig d/dq,) a (r, + Q^', qy) + j^^b+ (r, - (g, + Q^') , -g,)}, (f 9) 

i96+(r,-qi)/aT = [6+(r,-qi), = 
= - ^ {(1 + ^'A + - ^5 b+ (t, - (g, + Q^') , -qy) + -f^^a (t, q, + Q^', g^)}. 

The inclusion of umklapp processes leads here to a system of engaging equation of motion for magnon 
operators for different cells of reciprocal sublattice. 

Note that the first equation passes here into the second equation with the changing a {T,q±) ^ b (r, — q±), 
b —b and with use of the Hermitian conjugation. It is also not difBcult to find the equations for the second 
pair of Hermit conjugated operators a+ (r, qj^), 6(r, — q^). 

4 Diagonalization of the antiferromagnetic spin Hamiltonian 

Let us perform the unitary transformation (a variant of Bogolubov-Tyablikov transformation fRef.[17], sec. 13) 
to new linearly independent creation and annihilation operators (r, gj,, E±) and ^ (r, gj,, E±) for two types of 
magnon states, which propagate along direction a;-axis with dimensionless energy E± = E ±b, where E - is a 
continuous energy parameter, and wave vector component qy in the range of qy to qy + dqy. 

We will write out the following two of the four equations (two others are Hermitian conjugated): 

a(r, q±)- J [u {q^, E) aT,qy, E_) + v* (q^, E) ^+ iT,qy, E+)] dE, (20) 

b+ (T,-qx)= J [v {q^,E)ar,qy,E^)+u* (q±,i?)C+ {T,qy,E+)] dE. 
The new operators must obey the boselike commutation relations: 

(r, qy, E±) , e (r, q'y,E'±)] - <5 (g, - q'y) S [E - E') , 

[i{T,qy,E±),i{T,q'y,E'^)] = [e+(r,g„£;±),e+(r,g;,i?i)] = (21) 
= [e(r,g„i?±),€(T,g;,i?^)] = [e+(T,g„i?±),e+ (r,g;,i?^)] -O, 
and the following equations of motion 

id/dr^ (r,g^,i?±)-[e (r, g„ i?±) , /is (r)] = i?± ^ (T,g„i?±), (22) 

id/dT C+ (r, qy,E±)^ [t (r, g„ ^±) , /15 (r)] = ~E^ ^ {r, qy,E±). 

Let us next retain in the right sides of Eqs. p^ only members with Qx = = 0, that is in the com- 
mutation relations the umklapp processes between distinct cells of reciprocal sublattice will be neglected (N- 
approximation) . In this case in the right sides of Eqs. (fTO|) only terms with values Q" = are retained. 

With Eqs. (|2ip we will obtain in this approximation the first pair unitarity conditions for transformation 
coefficients: 

J [u iq^,E) u* (g^, qy, E) - v* (q^^E) v (g^, qy, E)] dE = S (q, - q'J , (23) 

J [v* {q^,E)u{q'„qy,E)~u{q^,E)v* {q'^, qy, E)] dE = 0. 

Note that in this case due to terms, which contain the derivative ig d/dqx, the equations of motion conserve 
the inhomogeneous external field dependence. The inclusion of umklapp processes would leads to corrections of 
higher-order for the parameter 5 <C 1. 
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Substituting Egs- ipU]) . ([2^ to Eqs. fT^ . and making under the integrals the coefficients ahead of operators 
^_ {r^Qy, E^) and {T,qy, E+) equal to zero, we obtain the system of four equations for determination of 
transformation coefficients u (q_L, E), v (q±, E) and their Hcrmitian conjugated values: 

[ig d/dq^ -[E-l- 6a)] u iq±,E) + 7q^w (q^, S) = 0, 

7q^M (q^, ^) - [^5 d/Oq., - {E + 1 + 6a)] v (q^, S) = 0, 

[igd/dq^ + iE+l + bA)] v*{q^,E)+j^^u*{ci^,E)^0, (24) 

7q^«* {q^,E) - [ig d/dq^ + {E - I - 6a)] u* (q^, S) = 0. 

Note that the first pair of Eqs.(IM| transforms in the second pair by replacing u {q±, E) v* (qx, —E), and 
V (qj_, E) — > u* {q±, ~E). In this case, the condition is violated. Hence, it follows that solutions with E <Q 
are not physical and should be eliminated. However, if the energy parameters E have the common positive sign 
in both equation pairs, they change from one to the other by complex conjugation and by following replacing 
of transformation coefficients: 

u{q^,E)^u*{q^,E), v {q^,E) ^ v* {q^,E) . (25) 



Therefore, we will deal next only with the-second pair of Eqs. ii^T:,, . 

Remaining next in the context of the N-approximation; we will multiply the first equation from Eas. (|19p to 
the left by (r, q^ + Qx^ qy) and the second equation to the right by 6 (t, — [q^ + Qx) , ~Qy)- Upon integration 
with respect q^ and summation over Qx, the expression for antiferromagnetic spin Hamiltonian (|15p takes the 
form 

hs = hs{T) = -Nil + bA)/2+ (26) 
+ E / {^a+{T,q^ + q^) da{T,q^)/dT^tdb+ (r, -q^) /9t 6 (r, - (q^ + Qx))}dqx. 

Substituting Eqs.d^Ol), (1221) to Eq.^, we now obtain 

hs (r) =-7V(1 + 6a)/2+ 

+ 51 / {[u*{'lx + Qx,qy,E)t{T,qy,E^) +v{qx + Qx,qy,E)^{T,qy,E+)]- 
Qx 

• [uiq^,E')EU{T,qy,E'_)~v*{q^,E')E'_^e{r,qy,E'_^)] - (27) 
- [v {q^,E) E^^ (r, qy,E^)~u* {q^,E) E+t (r, qy, E+)] ■ 
■ [v* {qx + Qx,qy, E') t {t, qy, E'_) + u {qx + Qx,qy, E') ^ (r, -B;)] }dq^dEdE\ 

where the transformation coefficients of type u* {qx + Qx, Qy, E) and v* [qx + Qx, qy, E') are derived by means 
of simple argument shift q^ — > + Qx,qy in coefficients u* {qj_,E), v* {q±,E). The spin Hamiltonian in the 
form ([27]) accounts, like Ea. (fT5)l . the magnon states correlations for different cells reciprocal sublattice. 

We will transform next the part of expression ([?7|) containing operator product £^ {T,qy, E+) £^ (^T,qy, E'_) 
through variables E' E permutation to the form: 

^ I [E'_v{qx + Qx,qy,E)u{q^,E')(,{T,qy,E+)i{T,qy,E'_)^ 
Qx '' 

-E^v{q^,E) u {qx + Qx,qy, E') ^ {T,qy, E^) ^T,qy,E'_^)]dq^dEdE' = (28) 
= J E'_dEdE'J2 J {qx + Qx,qy,E)u{q±,E')aT,qy,E+)^{T,qy,E'_)- 

-v{qj_,E') u{qx + Qx,qy,E)^{T,qy,EL)^{T,qy,E+)] dqj_. 

Let us require now the fulfillment of one more condition from the second pair of unitarity transformation 
conditions 

J [v{qx + Qx,qy,E) u{q±,E') -v{q^,E') u{qx + Qx,qy,E)]dq± ^ 0, (29) 
which reduces to zero Eq. . Instead of Eq. (|?7)) we will then obtain 

hs (r) =-iV(l + 6A)/2+ 
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-«* (g. + Q., i?') « (q±, i?) e (r, g^, i?-) C+ (t, qy. E'_)] + (30) 

+^+ [u* iq^,E) u (g, + Q,,qy,E') ^+ (t, g^,, S+) ^ (r, g^, S^) - 

-V* {q^.E) V (q, + Q^,qy,E') ^ (r, g,^, e+ (r, g^, } dq^dEdE' . 

Let us perform the magnon operators permutation in the second and the fourth rows of Eq. (|30p . Using the 
first commutation relation in Eq. (|2f p . taking into account that 

hm , 6 {qy - q'y) ^ (r, g^, E±) (t, qy. i?±)] = Nyl2TT 5 {E - E') (31) 

and also the second condition from the second pair of unitarity transformation conditions 

I [u* [q^ + Q, qy,E') u {q^,E) - v* (q, + Q.,qy,E') v (q^, E)] dq^ = 6 (E - E') , (32) 

we obtain for spin Hamiltonian (|27p diagonaHzed form: 

hs (r) --7V(l + 6A)/2- 
- Ny/ {2tt) J2 1 1^+^'* (q^' ^) « (9- + + (9- + Q^^%^ ^ (qi> ^)] ^^q^ci-s + (33) 

+ j [E+i+ (r, 9^, ^+) e (r, g„ S+) + {r, qy,E^) ^ (r, g^, S_)] 

The foUowing transformation of the second term in Eg. ([55]) will be given at the end of Section 5. The last 
term here describes magnon excitation energy. It is consistent with the above-considered equations of motion 
for magnon operators (1^^ . Because the quadratic form describes a stable state, close to ground state, it 
should be positively definite. That is, the energies of both types magnons should be E± = E ± b > 0. 

5 Asymptotic solution of equation for the unitary transformation 
coefficients 



Let us go now to calculations of the transformation coefficients. We will use the second pair of equations ll24|) . 
in wich the value v* {q±, E) will be eliminated. Because the variable qy in Eqs. ([24|) is considered as a constant 
parameter it will be not explicitly indicated in arguments and the new variable will be used 



c 



j Ici^dq, (34) 





Then we will switch from the notations u* {q± , E) and to u* {(, E) and . The partial derivatives 
with respect to ^ will be replaced by the ordinary derivatives. As a result, we obtain the following differential 
equation for the transformation coefficient u* (C, E) 

g^d^u* (C, E) /de + P (C, 9) gdu* (C, E) /d^ + r (C, 5) u* (C, E) = 0, (35) 

where 

p(C,5)=Po(C) = -2»i?/7C: (36) 

r (C, 5) - ro (C) + 5 ri (C) - [(1 + 6a)' - E^ - 7^] + {1 + Ba - E) -d (I/7,) /^C- 

Solution of Eq. psp will be presented in the form of an asymptotic expansion in terms of small parameter 
5< 1 (Ref.[18], sec. 7.1.6): 

00 

u*{C,E) = Y,9"iAn{C)expei+Bn{C)exp92), d9,/dC ^ {() /g, (37) 

n=0 
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where the variables Oi, C are assumed independent from one another and (C) will be defined below. 
The derivatives in Eq. (l35|) are transformed as follows(A^ = dXi/dC,): 

g d/dC = -Xid/dOi - X2d/d02 + gd/dC, 

d^/dC^ = Xjd^/del + 2X1X28^ I deide2 + Xld^/Oel - (38) 

-2gXid^ / dOidC, - 2gX2d^/de2dC - gX^d/dOi - gX'2d/d02 + g^d^/dC^. 

The substitution of expression ([57]) in Eq. ([55|) and making the expressions ahead g"exp0i_2 equal to zero 
gives two equations 

{Xl + Aipo (C) + r (C, g)) A„ - g (2Ai + p (C, g)) K ~ gX[A, + g^A"„ = 0, (39) 

{Xl + X2P0 (C) + r (C, g)) B^-g {2X2 + p (C, 9))B'^^ gX'^B„ + g^B'',^ = 0. 

By assembling then in Eq. pop the coefficients near the same powers of we will derive the equations for 
determination of values An (C), S„ (C). Coefficients ahead g'^Ao g°i?o (C): (^O: Bq ^ 0) give the equations 
for determination of A^ = Ai 2: 

A2+po(C)A^+ro(C) =0, 1 = 1,2, (40) 

from where it follows 



Ai,2 = */7c(-^^±^^(C)), E{C)= ^{l + bAf-jl (41) 

The zero-order quantities Aq (C), Bq (C) are determined by differential equations obtained by making the 
coefficients ahead of the first power factor g in Eq. (f39|) equal to zero (we will restrict next only by the first 
order): 

(2Ai+po)Ao + (A'i+ri)Ao-0, (42) 
{2X2 + po) B'o + {X'2 + ri) Bo = 0. 
Solutions of equations ((1^ have the form 



Aq (C) exp 9i=A exp h (C) , (C) exp 62 = B exp h (C) , (43) 

where 

/i,2 (C) = - / l^Y^^dC -- j X,,2d( , (44) 
J 2Ai,2+po .9 J 


2Ai^2+Po-±2zS(C)/7c, (45) 
^ ""'-dC = -1/2 d [In (E (C) /7c) T ( C) /tcI (46) 



2Ai,2 +Po 



and constants of integration A = Aq (0) = A {E, qy) , B = Bq (0) = B (E, qy) depend here on E and qy, as on 
parameters. 

Due to i5 > the quantity Ai (the upper sign in formulas (|41]l . (|45p . ((46|) ) may have any small values. In 
the limit of 5 ^ we have Ai 0, that is in the case of homogeneous field we obtain the known result or two 
branches of magnon spectrum: 



E±=E±b^E± (q^, b)=EiO± b= ^{1 + bxY - 7c ± & > 0, (47) 

whereas the going to the homogeneous case in expression for A2 leads formally to nonphysical negative value of 
parameter E < Q and conditions ([23|) . ([32|) are violated. We will choose out of two approximate solutions the 
solution u* (C, E) sa Aq [Q) exp 9i = A exp Ii (C) that corresponds to low branch of magnon spectrum with gap 
= S (C) - 6 > 0. 



This gap disappears at the condition E (qj_ = 0) = y (1 + ^a)^ — 1 = fee ~ a/ 26a + b\ = b, which is 
consistent with known condition for orientational phase transition of spin- flop type in antiferromagnet at critical 
field b = be fRef. [19j). The equilibrium state of two neighboring spins, belonging to different sublattice with 
zero total spin polarization passes into state with nonzero total spin polarization (Fig. |4|) . 



Taking into account now that E {() = E {q±) = y 1 -I- 5^ — 7^^, E { ( = 0) = E (qy), we will write the 
exponent Ii (C) from Eq. ljlSjl in the form 
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Figure 4: The scheme of spin-flop phase transition in easy axis antiferromagnet at critical field b = be- For 
FeCOg: cos 6* = bc/2 « 0.22, 29 « 155°. 



-l/21n||^^ + l/2 [i?(q^)/7q,-i?(g,)/7,J-^|i?fe- /i?(q±) dq^ 



(48) 



To obtain the approximate expression for Ii (q^ , E) , we will perform for variables the extension trans- 
formation of type qx/g = V (Ref.[l8J, sec. 4.1). In the limit, g ^ the "extended" variable takes values 
< \r]\ < oo and asymptotic expression for Eq. (|i5|l takes the form 



Eri- j E{gi], qy) dr] 



(49) 



Turning next to the earlier "no extended" variables, we will obtain the following approximate expression 



u* {q±,E) = Aexp 



- ( Eqx - / E{q±_) dq^ 







0(ff) 



(50) 



Using next the recent equation in system ([M]) for the second coefficient v* {rj, qy, E) we will find the following 
asymptotic expression 



V* iq^,E) = - — 

7q_L 



1 

7q, 



ig d/dq, + ( ^l + bl - E 



u* (qx,S) 



(51) 



-^ll + ^c + E 



Aexp 



^ f 

- I Eqx - J E{q±) dq^ 



0{g). 



The calculations of coefficient A are placed in Appendix Al. We have 



1 I ; 

47r5 I E 



(52) 



Finally, for the unitarity transformation coefficients we will obtain the following expressions: 



u* {q^,E) 



Zing V E 



exp 



cxp 



Eqx - J E{q± ) dq^ 



Eqx- J E (q^) dqx 



Oig), 



0{g). 



(53) 
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The integration over variable E in the second term of the right side of Eq. ([S5|) with the help of Eq.(Al.4) 
now gives 

J [E+iv*{q^,E)viq^ + Q,,,qy,E)) + E_ {v{q^,E)v* {q^ + Q,,qy,E))]dE ^ 



1 VTTb\ 



E 



Qx_' 



-{E-E{q^))Q, 
9 



d{E-E{q^)) + 0{g) 



^i + bl - E (q^)^ S (g,) +0{g)= {j I + bl - E (q^)^ N.J2^ + O {g) , 



where we accounted that 



5{Qx)= hm 5{q^- q^ - Q^) ^ N^/2ii5Q^fi^ N^/2tt . 



(54) 



(55) 



As a result, the expression for the diagonalized Hamiltonian ([33|l takes the following approximate asymptotic 
form: 

hs = -^l + bl iN/2 + N^) + N^l {2TTf J E{q^)dq±+ (56) 

+ J [E+t (r, qy,E+) ^ (r, qy, E+) + E^^+ (r, qy, E^) ^ (r, qy, E^)] dEdqy + O (g) . 

The first two terms in the right part of Eq. (|56p describe the ground state energy of antiferromagnet with 
regard to zero-point oscillations. The third term corresponds to the energy of zero oscillations inhomogeneous 
two-dimensional spin system. 



6 Indirect interspin interaction close to antiferromagnetic spin-flop 
phase transition 

Let us set off two substituted atoms with nuclear spins at sites k and I from others of substituted atoms in 
one-dimensional chain. They are arranged in the plane of plate Zk = zi = and separated along x-axis by 
xik — xi — Xk ^ I — k distance. For simphcity, we assume that these atoms belong only to sublattice A (Fig. [3]). 

Due to the interaction with antiferromagnet spin system the nuclear spin of fc-th atom excites a virtual spin 
wave, which is propagated along the plate and absorbed by nuclear spin of other atom /. As a result, the indirect 
interaction between them is produced. For the case of homogeneous external field, the indirect interaction was 
considered early in Ref.[8-10]. The contribution of magnetic dipole nuclear interaction with electron spins of 
neighboring atoms was considered in Ref . [11] . 

We will here account the inhomogeneity of external field with weak gradient along one-dimensional nuclear 
spin chain, which represents one-dimensional quantum register. Hamiltonian of interaction of k-ih nuclear spin 
with external magnetic field and of hyperfine interaction with electron spin of own atom will be described by 
the following expression: 

His [k] / (^wb) = his [k) = ~m (k) /a. (fc) + a Ia (k) Sa (k) = (57) 

= - (u,i (fc) - a^A. (fc)) /a. (fc) + a/2 [/+ (k) (k) + (fc) S+ (fc)] , 

where Iaz, I a — ^Ax i ilAy are components of nuclear spin operators, lui (k) — {^i /js) {b + gk) is resonance 
nuclear frequency in local field bk = b + gk, ji/^s ^ 10~^, a — A/ue ^ 10^^ is isotropic dimensionless constant 
of hyperfine interaction, A/2tt ^ 100 MHz. 

Let us go now in Eq. (j57p to interaction representation. Keeping in mind that antiferromagnet is considered 
in the very low temperatures (close to ground sate) we will retain only linear dependence from electron spin 
operators {Saz ~ 1/2): 

his (t, k) = ~ {lui (fc) - a/2) Iaz (k) + Ahis (r, k) , (58) 

where Hamiltonian of nuclear and electron spins transverse hyperfine interaction of the fc-th atom, considered 
as perturbation, has the form 

Ahis (r, k) = exp (z/i^r) Ahis (fc) exp {-ihsr) = a/2 [l+ (k) (r, k) + (fc) (r, k)] . (59) 
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Taking into account Eqs.(f51l.(fTT |) . (P^ . for the considered quasi- two-dimension antiferromagnetic structure 
we will have (in the following, index A will be omitted) 

S- {r,k) = {S+ (T,fc))+ = 

1 X - 1 f 

a+ (r) exp (iq^k) -> -— — / a+ (r, qj_) exp (iq^k) dq^ = (60) 



q 



= / ["* (^-L' -E^) exp {iE^t) C+ {qy, E^) + V (q±, S) exp {-iE+r ) ( {qy,E+)] exp (iq^k) dEd(i±. 

Let us assume that hyperfine interactions of nuclear and electron spins of own atoms k and I are switched on 
adiabatically slowly at r' = — oo, when the unperturbed density matrix for two nuclear spins in antiferromagnet 
is represented in the form of direct matrix product p{—oo) = pi {—oo) ps {—oo), and acts to moment t' — t. 

We will write the equation for electron density matrix ps (t') in interaction representation 

tdps (r') /dr' = [{Ahis {r' , k) + Ahis (r', 0) exp {-s\t - t'\) , ps (r')] , (61) 

where s is a small parameter, which characterizes here the rate of interaction switching (s <C 6c ~ ^) in the 
adiabatic development processes from initial density matrix ps (— oo) to matrix ps (r). It assumed the value of 
parameter s is order of the magnon damping rate or width of antiferromagnetic resonance line, which for ideal 
crystal at low temperatures may be less than or of the order of several oersteds (~ 10"^ T) (Ref.pJ]). 
It follows that in the first order of perturbation theory from Ea. (|6ip we will have 

T 

PS {t) ^ PS {-(x) ~ i I [{Ahis{T',k) + Ahis{T\l)),ps{-^)]e^v{-s{T-T'))dT\ (62) 



where in the limit of low temperatures only matrix element of non-perturbed density matrix ps (— oo) for ground 
pure state |0)(0| is nonzero, that corresponds to the absence of antiferromagnet magnon modes: 

m'^{qy.E±)a<ly,E±)\Q)=Q. (63) 

Using now Eq.(|62p and taking into account that partial trace over electron spin states 
ti's Ps (~oo) {^his (r, k) + qAhzs {t, 0) — 0, we will obtain the second order correction to Hamiltonian of two 
nuclear spins as the mean value of their hyperfine interaction over electron ground state: 

hii {k, I) = trs PS {t) {Ahis {t, k) + Ahis (r, 0) ~ 

r 

^-i J2 J [Ahis{T\j), ps{-^)] Ahis{T,j')exp{s{T' -T))dT' = (64) 





[ {0\[Ahls{T,j),Ahls{T + T',j')]\0)e^p{.ST')dT'. 



j,j'=k,l_ 

Taking into account that (0|S'='= (r, j) , {t + T',j') |0) = 0, we wiU write 



hii (fc, I) « —ia^ ' 



Jj'=k,l_^ 

After cumbersome rearrangement of Eq. (165p (see Appendix A2), we will obtain the expression for indirect 
interaction between two separated nuclear spins, belonging to common sublattice: 

U{k,l)^ (66) 

u* {q^,qy,E)u{q'^,qy,E) v{q^,qy,E)v* {q'^,qy,E)' 



2(27r)^ 



Re 



E^ -I- is E+ — is 



exp [i {q^k — q'^l)] dEdq±dq'^, 



Let us retain now in Ea. (|66p only major members, containing in dominator E^ = E — b (the energy of the 
low magnon mode) and change the variables of integration as follows: 

<lx,q'a: {qx - qx) = q-, {qx + qx)/2^q+, dq^dq'^ = dq-dq+, 
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(g^fc - q'J) ^q-{l + k) /2 - q+ - k) (67) 

Once again, we will use the transition in the exponent to "extended" variable w = q-/ g and go to limit <? — > 
in the preexponential factors that is we will take q- — q^ — q'^ = gw. This allows to perform the following 
changing g+, q'^ ^ q^, and dq^dq^ gdwdqx- As a result, we will have in place of the expression 

U{k,l) = 



2 r l„, /'^ |2 



= -^gRel l^(q..^)Pexp [^jE-Ejc,,) +.g(^ + fc)/2) - ^qAl ~ k)] 

Note here, that the value of denominator ~ E — h in (|68p may be small but larger than some minimum 
value, required for implementation of permutation theory condition, which is determined by ignored weak 
nonlinear effects. 

By extending limits of integrating over w in Eq. (p5|l to ±00, we will obtain 
00 

j {exp [i{E-E (qj.) +g{l + k) /2) w]}dw « 2t:5 [E - E (q^) + gk + g {I - k) /2) . (69) 
—00 

By integrating then over E and using Eq. (|53p , as a first approximation of asymptotic expansion in terms of 
(7^1 we will have 



y^l + bl + E{q^) ^ COS {q^ {I - k)) 
4(27r)^"^y E{q±) E {qj_) - \i+k)/2 + ii 



U (k, I) = ^Re / v-^^cj^vMx. . + Oig), (70) 



where 6(/+fc)/2 = b + g {I + k) /2 = bk + g {I — k) /2 is the local external field in the middle point {I + k) /2 of 
l-th. and A;-th atoms. It is further conveniently to choose the field value bk = b + gk aX k-th site as the origin of 
inhomogeneous field reading. 

In the following, we will simplify the expression for magnon spectrum 



E{ci^) = .^l + bl-^l^. (71) 

For this we will replace the initial q_L dependence of for all components of wave vector values < Ifel, < 
TT by parabolic isotropic dependence. That is we perform in our model the following replacing 

7^^ = [1/3 (cos (fe/2) + cos {qy/2) + l)f -.^ 1 - qi/12, (72) 

wich coincide with the expression received by simple expansion for small values q'^ of initial spectrum: 

7^^ « 1 - qi/12 + . . . . (73) 

Because the small values of q^ are the most essential for the following, we will use the form ([7S|) . extrapolating 
its in the expression under integral over all values < qj_ < tt. 

Using this simplification in the expression ()71|) . we will go to the new variable of integrating C + = 
E {ci±) ~ \/bc + qi/12 and dq± = q±dq±d(p = 12 + be) d£^dip. As a result, we will obtain (the symbol 
0{g) will be then omitted) the expression, generalizing the known Nakamura's expression for antiferromagnet 
in homogeneous field [7]: 

^U{k,l) = V{Abk,l-k)= (74) 



^bl+.yi2-bc (^TT&i: + 6c + e) (e + Afefe -9(1- k) /2) Jq (^^ 12 ({^ + be f -bl) [l - k)^ d^ 



{i + Abk~g{l~ k) I2Y + s2 

u 

where the integral over ip is expressed in terms of zero-order Bessel function of the first order (Ref.|20j. Eq. 3.715- 
18.): 

2-7r 



j cos [qj^ (/ — k) cos if] dip — 2tt Jq (q^ {I — k)) (75) 



and it is denoted 

,2 



245c {bc-b-g{l + k) /2) = 246c (A6fe ~ g {I ^ k) /2) = M(i+fe)/2 , A6fc = 6c - 6 - gk. (76) 
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For M(;+fc)/2 ranges of interqubit distance are 

1 < l-k < 2Abk/g. (77) 
For M(i+fe)/2 = -'^fi+fc)/2 < the ranges are 

2Abk/g<l-k (78) 

The critical value of the middle point position (Z + fc) /2 = be or A6fc = g {I — k) /2, wherein the parameter 
l^\i+h)/2 changes sign and the dependence of indirect interaction of nuclear spins from interspin distance is quali- 
tatively changed, we will denoted as "turning point" , by analogy to problem on the barrier reflection of quantum 
mechanical particle in WKB-approximation (Ref.|2I], sec. 46). Note, that further it will be conveniently to tune 
through parameter IShk — be — b — gk the state of nuclear spin k in quantum register and to use spins I in the 
turning points, separated by distances I — k ^ Luk = 2Abk/g, where L is period of chain, Uk = 1, 2, 3, . . .. To 
perform the two-qubit operations in quantum register it should switched on the interaction between considered 
nuclear spins. The way to do this is through the tuning of the spin-qubit state to the turning point state. This 
requires the fast as compared to operation and decoherence rate variation of external field. 

Let us account now that major contribution to the integral in Eq. (j74p is given by values ^ ^ 1, s <C Abk 
and tt'^ ^ /^(i+fc)/2 whereby the denominator is small. Let us retain next under integrals only quadratic on 
qj_ terms and assuming then the rapid convergence of the integral on upper limit, we will replace tt — > cx). As 
a result we will obtain for indirect interaction the approximate expression (Ref. 20 , Eqs. 6.532-4, 9.561.5): 

^U{k,l)^V{AbkJ-k)^ 
ri" , , \ /■ Jo(q±(Z-fc)) q^dq± f Jp {q± {I - k)) q±dq^ 



•^_L + M(;+fe)/2 

^1 + bl + bc^ Ko{fiii+k)/2\l~k\) , 

where Kq (x) is zero order Macdonald function. 

Its asymptotic value for t: {I — k) > /i(;+fc)/2|^ — A;| ^ 1 is: 

V (Abk,l - fc) « 2 [^Jl + bl + 6c) y ^^^^^J^|^_^ |exp [- (M(i+fe)/2l^ - k\)] .. (80) 

As is seen from Eq. ()80p . the parameter (/^(;+fc)/2) ^ presents the effective radius of indirect interaction 
between nuclear spins. In the case of homogeneous field it rapidly increases close to value for bulk spin-flop 
phase transition, when be b. In the case of inhomogeneous field effective radius increases at condition 
l^\i+k)/2 ~^ t'^^t when only local field in the middle point for the considered two spins is close to the value 
for bulk spin-flop phase transition b + g {I + k) /2 ^ be or Abk + g {I — k) /2 — > 0. At the same time, relation 
be > b takes place and homogeneous phase transition does not occur. Note that if aforementioned condition is 
not fulfilled, the radius of indirect interaction remains of the order of lattice parameter. 

To receive the real value V {Abk,l — k) for the case M(;+fe)/2 — ~'^fi+k)/2 ^' the lower limit of the first 
integral in Eq. ([79| should be replaced by W(;+fc)/2 > 0. As a result, we will obtain (Ref. [20], Eqs. 3.753-4, 
8.405-1, 8.407-2): 

V {Abk, l-k) ^2 I^^Jl + bl + be^ ReXo {-i'^[i+k)/2 {I - k)) = (81) 

= -2(^\/l + 6^ + 6c) ^No{iyii+k)/2{l~k)), 

where Nq (x) is zero order Neumann function (or zero order Weber function Yq (x)). Using its asymptotic value 
for V(i+k)/2 il-k) = ^2Ahe {g {I -k) 12- Abk) (l - k) ^ ^ 2Abcg {I - k)^ /2 > 1 we wifl have 

V{Abk,l-k)^2(^^l + bl + be^ y 2^(;+ J2|^ - fc| ^''^'+''^/^'^ ~ ^' " ^'^^ ■ ^^^^ 

In this case, dependence of effective interaction from distance between spins takes oscillating character with 
the quasi-period, which is decreased with increasing of the distance I — k. As this take place, the interaction 



goes to zero at the sequence of distance y 24beg {I — k) /2 « (n 4- 1/4) tt 
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Figure 5: Dimnsionlcss indirect interaction V (A6/5, 1 — k) as a function of interqubit distance I — k for "exact" 
expression (|74|) (A) and for approximate analytical expressions ([79|) . ((8T|) (B). The used numerical values of 
parameters are: h^, = 1/4, .g = 2 • 10"^, L = 100, s = 10"^ and A6fc = 6c - 6fe = 1 • lO^^, 3 • IQ-^. Note, 
that transition to value g ^ 2 ■ 10^^ corresponds here to the following modification of turning points positions 
(vertical lines): Luk = 2/S.hk/g — > Lnk ■ 2 • /g. 
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As is seen in Fig.[5l in the case that value of damping s — 10^^ ^ A&fc, the approximate distance dependence 
of indirect nuclear spin interaction (in the form of Eqs. ([75|) . (|8ip ) outside of the turning points is in good 
agreement with exact relation fEq. ([7i)) ). 

In the limiting cases for fJ-(i+k)/2 (^ ~ ^) ^ 1? '^{i+k)/2 {I — k) ^ 1 we will have, respectively, the following 
logarithmic growing with decreasing \l — k\ asymptotic dependences of indirect interaction (Ref.[20J, Eqs. 8.447-3, 
8.444-1) (C = 0.577. ... is Euler constant): 

V {Abk,l ~k)^2(Jl + bl + bc] In (2/ Ua+k)/2 {I - k))) , (83) 



k) « 2 (^^l + bl + bc^ In (2/ (a.(;+,)/2 (/ -k))), 
k) « 2 l^^l + bl + bc^ [In (2/ (/ -k)))- C] 



V {Abk, l^k)^2l^^l + bl. + bcj [In (2/ {iyii+k)/2 {I - fc))) - C\ . (84) 

Note however, that using the second order of perturbation theory for calculating the indirect interaction 
V (A&fc, I — k) immediately close to the turning point is not sufficient and the account of magnon damping and 
relaxation effects are required. Therefore, minimal values 'y(i+fc)/2 and fJ'(i+k)/2 should be bonded by same small 
nonzero values. This values will be bound also the range of maximal distances wherein the expressions (|79l) and 
(|8ip are applicable and, respectively, the size of quantum register is estimated. 

We will obtain the corresponding expression for U {k, k) assuming / = k and s ^ Abk in Eq. ([7i|) : 



y/b%+r^-/12-bc 



i-2Uik,k) = ViAbk,0)^ J ^ e + A^. ' ^ ^''^ 



T7^, + bc- Abk) m + ^ -^ + ^^7^2- be. 

It has relatively weak logarithmic dependence on the position of fc-th nuclear spin. 

Let us estimate, finally, the value of direct dipole interaction between nuclear spins-qubits in quantum 
register. Hamiltonian of transverse part of dipole interaction for nuclear spins k and I has the form 



where dimensionless ratio 



in which 



is the dipole field produced by nucleus, which is located at a distance of the order of sublattice period Ux ~ 
1 nm, = OAtt T^cm^/J is vacuum magnetic permeability. This field is of the order of 10"^ T, the ratio 
D - 10-1" < a2 - 10-^ 

The ratio of dipole interaction to the indirect interaction (in the form of Eq. (|80p is: 

D _ V2^D tii+k)/2 ^gg^ 



[/(fc,0 |fc-;|3 3„2(yi^ + 6^j i?5/2exp(-i?)' 

where R = A*(i+fc)/2|^ ~ k\ > 0. 

Maximum of R^^^ exp {—R) obtained at i? = 2, 5 is 0,78. For the same value of R and for distance |fc — ^| > 
100, parameter /X(i+fc)/2 corresponds to value < H{i+k)/2 <C 1 or < 2'ibc {be — b — g {I + k) /2) <C 1 . In the 
case, when the value of local field in middle point is close to the critical field of spin-flop phase transition, the 
indirect interaction becomes very large and dipole interaction for considered nuclear spins would be negligibly 
small in comparison with indirect interaction (Fig.[6l). It follows that a possibility of considerably increasing of 
interaction arises for a definite pair of removed qubits as compared with interaction between all another pairs. 
Furthermore, if the local field in middle point is more than the critical field or 2Abc {be — b — g (l + k) /2) < 0, 
the indirect interaction takes oscillating character (see Fig.jHl). 

Let us assume for example, that register has 1000 nuclear spins-qubits that is the total distance |/ — fc| 
in register with period 100 qubits is of the order of / — fc = 10^ and value of M(;+fe)/2J wherein the indirect 
interaction is essential, corresponds to turning point or Abk ^ be — bk ^ g {I — k) /2 ^ \. 
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Figure 6: The indirect interaction V (Afefe, I — k) is large for qubit pairs I, k marked by the local field close to 
middle point 6(fe+;)/2 ^ be or for 1 ^ l^\i+k)/2 > At the same time, this interaction is far less for other qubit 
pairs with b(k+i)/2 < be- 

7 Indirect interspin interaction close to antiferromagnetic resonance 
(AFR) 

The homogeneous antiferromagnetic resonance (AFR) phenomenon refers to the excitation in antiferromagmet 
by microwave field's homogeneous mode with q = 0. Resonance in thin plate will be well resoluted when the 
energy difference of principal homogeneous transverse mode with qz = and of the next transverse modes with 
Qz = 27rOz/(i, will fit the condition 

^Jbl + {2Traz/df / 12 - be > Aujo/lJe, (90) 

where Awq is AFR width. It follows that be > (27ra^/c?)'^ /246c' > Acjq/'^b that is thickness d of large plate 
and correspondingly the site numbers of sublattice N^, should be determined by inequalities 



^ <Nz^d/az<nJ-^, (91) 



6bc ^ ^ V 66(7 Awo' 



'Jc 

Consequently, for example, if Aloq/uje ^ W^^bc the thickness should be in the range of micrometers. Note, 
that the condition permitting to restrict by only zero mode in the case of Sec. 2 differs from the condition (|9ip 
by the change of the AFR width to the quantum states decoherence rate 1/Td. 

Let us assume that the antiferromagnet plate interact with the homogeneous transverse microwave field 
left-hand circularly polarized relatively to external field (dimensionless amplitude <C 1 and frequency lu): 

bx (t) + iby (r) = b± exp {iujr) . (92) 

The interaction Hamiltonian in interaction representation, accounting that directions of sublattice spin 
vectors have opposite sign (Ref.[17j, sec. 4.2), has the form: 

N 6 

AHs (r) / {hcjE) = Ahs (r) = b^ exp ( hot) ^ {S+ (r, r,) - 1/6 ^ S+ (r, (r, + r^))} + H.c. « 

i=l S 

«6j^exp (iwr) ^/n~[ V,Qx^9„,o (flq^ (t) -7q_L^iqi W) +H.C. ^ (93) 



where 



^ J F{E)exp{-iA{E,uj)T)({0,E^) dE +H.c., 
F(£;) = 2^&i^[7.(Q,,0,^)-7Q,z;(g,,0,^)] = 



= 2TThi_ J2 (0, E) - 1/6 (exp (zQ,/2) + exp (-iQ,/2) + A) v (0, E)] exp (-z {E - be) Qx/g) ~ (94) 

Qa, 



U^S{E-bc) + 0{g) 
V V 1 + 
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Here the relations (O, (fTTj) . ([^O]) are used and it is accounted only lower branch of magnon modes with 
energy To receive the final result in the Eg. ([94|) . we have used Eq. (fT3|) for 7q^, Eq.(Al.4) for sum over 
Qx — 27rm and accounted also relation g/2 ^ be < I- It follows from Eq. ([93|) that magnon modes can be excite 
only with Qy = 0. The difference A {E, ut) — — uj~E — b — ut plays here the role of energy gap for lower 
branch of magnon spectrum in the frame of rotating with microwave field reference. In the absence of gradient 
of external field, this gap is the frequency detuning of microwave field relatively to frequency of homogeneous 
AFR: ajQ = be — b. It follows that AFR may be considered as a spin-fiop phase transition in rotating frame of 
reference at critical field be = beff {ijj) = + At the same time, the equilibrium phase transition at this 
field does not take place if be > b. 

Let us use now the so called resonance approximation fRef.[21j. sec. 40), whereby the frequency detuning 
is assumed to be small |A {E,uj) \ <C OJ and accounted only two magnon states. One state corresponds to the 
homogeneous non-perturbed ground state |0) in the absence of magnons, when (Ol^"*" [q'y, E'^ {qy, E)\ 0) = . 

SE Sqy 

The second state J J \qy, E)dEdqy corresponds to the excited magnon packet state in rotating frame of 



reference. It is characterized by the narrow energy interval SE <^ be near E = E {q± — 0) — be and by 
the minimal discrete interval Sqy — 2TT/Ny for wave vector component qy near qy — 0. The ground and magnon 
pocket states will be normalized as follows: 

SE Sqy 

(0|0) = 1 and J J {qyE\qyE)dEdqy = 1 (95) 



Let us represent by microwave field perturbed state vector |r) as superposition of this two steady states: 

SE Sqy 

\t) = |0) Co {t) + j J \qy, E) ci {E, qy, t) dEdqy. (96) 



We will next consider only two non-diagonal matrix elements for perturbation operator ([93]) , which have the 
forms: 

(0\Ahs (r) \qy,E) = exp (-z (be - b - uj) r) J F {E') (0|C {qy, E'S) \qy, E)dE', (97) 

{qy,E\Ah+ (r) |0) = exp {i {be - b - u) r) J F* {E') {qy,E\e {qy,E'_) \0)dE', 

where the frequency detuning value in this interval will be taken constant: | A {E, uj) \ ^ be — b ~ uj. 

By using the Schrodinger equation for interaction Hamiltonian (|93p (magnon relaxation processes for sim- 
plicity are here neglected) 

id\T)/dT ^ Ahs{T)\T) (98) 

and accounting the conditions (I95p . we will obtain the system of equations for determination of amplitudes 
Co (r) and ci {E,qy,T) w ci {be,T): 

SE Sqy 

idco (r) /dr = b^exp [-i {be-b-Lu)T] J F {E') dE' J J {0\^ {qy, E'_) ky, E)ci {E, qy,T) dEdqy, (99) 



idci{E,qy,T)/dT = b^exp[t{be~b^Lu)T] j F* {E") dE"{qy,E\i+ {qy,E'L) |0).co (r) , 
which should be solved with the normalizing condition 

SE Sqy 

\co{t)\^ + J J {qy,E\qy,E)\ci{E,qy,T)\''dEdqy^\co{T)\'' + \ci{be,T)f^l (100) 



and initial condition |co (0) P = 1. 

By eliminating next the amplitude ci {E,t) from system ()99p . we will obtain for determination of cq (t) the 
following equation 

d^co (r) /dT^ +i{be-b-uj) dco (t) /dr + \Bfco (r) = 0, (101) 
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wherein 

SE Sqy 



\B\^= j F{E')F*{E")dE'dE" j J m{qy, E'_) \qy, E) {qy, E\t {qy, El) \0) dEdqy = (102) 



= J FiE')F* {E")dE'dE" 5qym{0,E'_)e {p.E'!_)\Q). 
Let us use now expression (A2.4) 

(0| [e {qy,E'_) , t {qy,El)] I 0) = (0|e {qy,E'_) ^ K I «) = {iv - <l'y) ^ - ^") 



and go according to Eqs.(6), (11) from continuous to discrete values qy, namely, [qy, E-) \/Ny/2Tr £^q^ (E-) 
and 5 {qy — q'y) Ny/2TT Sq^^g'^ (compare with Ea. ((?T|) ). Taking then qy = and using Eq.(A1.6), we will 
obtain for the matrix element product in Eq. p02p the expression: 

hm m {0,E'_) e {0,E'l) |0) ^ Ny /27r lim (Ol^o {E'_) ^+ {El) |0) = (103) 

hi — —DC h, — "fh. —be 

= Nyjl^ lim b{E'- E") = Ny/ (27rg) . 

E'^E"=bc 

As a result, by using then Eqs. (l94|) . p02p . we obtain 

\B\' = ^^ + 0{g). (104) 

Consequently, the weak inhomogeneity of external field {g ^ 0) practically does not influence the value of 
parameter i? 7^ 0, which goes to zero in the absence of anisotropy (6c = 0). 

By using Eqs. (|100|) . ()101|) . we will now obtain the following expressions for coefficients cq (t) and c\ {be, t): 

/X / , X / rt{uj)T i {be — b — uj) . fl{Lu)T 
Co (r) = exp (-Z {be - b - uj) t/2) f cos ^-^—^ ^ sin 

ci {be, t) = i exp {i [be-b- u) t/2) sin (105) 

The probability that the antiferromagnet is in ground state at the time r will be determined by expression 

that is antiferromagnet ground state at homogeneous AFR is modulated with Rabi frequency 



n {cj) = ^J{be-b~Ljf +4\B\^. (107) 

At exact AFR resonance {be — b = u), the probability of being in ground state 

|co (r) |2 = 1 - sin^ {\B\t) = cos^ {\B\t) (108) 

varies from unity to zero. The elements of density matrix, corresponding to ground state in AFR condition, are 
transformed to the form 

|0)(0| - |co (r) no)(0| = f 1 - ^^'^1' (1 - cosn{u;)r)] |0)(0|. (109) 

\ {be - b - ivy + 4\B\^ J 

To obtain the additional part of indirect interaction AV {Abk, I — k) caused by this modulation (the second 
term in brackets of Eq. (|109p we will change in denominator of Eq. the energy of low branch E — b into the 

more low- lying energy E — b—V, (w) and multiply the matrix element of ground state by factor ^^^_^_^ja^^|g|2 ■ 
As a result, we will have 

AV{Abk,l-k)K (110) 



2\BmJTTyh + be 



{be - b - ujf + ilB]^ 



[V {Abk, l~k)~l/2V {Abk -n{Lo),l~k)-l/2V {Abk + f7 (w) , ^ - fc)] , 
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The first factor in Eq. (lllOp has sharp maximum at frequency of homogeneous AFR. The role of gradient of 
magnetic field appears mainly in dependence of indirect interaction on the distance between nuclear spins. 
The positions of additional tuning points are defined here by following expression 

Abk -g{l-k)/2 = ±^{bc-b~ uf + 4|B|2. (Ill) 

As is seen from Eq. ljllOp . amplitude <C 1 and frequency lo of microwave field are the additional parameters, 
controlling the indirect interqubit interaction. At constant external field and microwave power the new tuning 
points is achieved at microwave frequencies 

uj±^bc-b±^ (Abk ~g{l~ k) /2f - 4|B|2. (112) 

Note, that qubit tuning to one of this turning points states may be performed, as distinct from relatively 
slow external field variation; by fast variation of microwave field frequency that has also of direct interest. 



8 Nonadiabatic decoherence and longitudinal relaxation of one qubit 
states in antiferromagnet-based NMR quantum register 

Let us consider the decoherence and relaxation processes of quantum state for single nuclear spin at position 
k on axes x of sublattice A are caused by its interaction with virtual magnon excitation in antiferromagnet. 
These processes will be described by transverse and longitudinal relative to external field (quantization axis) 
components of Bloch vector (Ref.|llj) 

P- {k, t) = {P+ (fc, r))* = (fc, r) - iPy (fc, r) - 2 tr {l' p (fc, r)) = 2 tr/ [l' pi (fc, r)) 

F,(fc,T) =2tr7(/,p/(fc,r)), (113) 
tr/42 = l/2, tr//+/" = l, 
where non-steady reduced nuclear spin density matrix in antiferromagnet is determined by 

PI (fc, r) = trs p {k, t) = 1/2 [1 + 2P, (fc, r) + P- (fc, r) /+ + P+ (fc, r) /-] . (114) 

The interaction of fc-th nuclear spin with external magnetic field and magnon excitations will be described 
here by Hamiltonian with isotropic hyperfine interaction constant of the form 

h^hs- {lui (k) - a (1/2 - a+ak)) hz + a/2 {l+ S + S+) = (115) 
^hs- {uji (fc) - a (1/2 - {Q\a+ak\Q))) hz (k) + Ahis (k) , 

where the average value of flopped electron spin number in ground state of antiferromagnet on fc-th site known 
as "spin contraction" has the form 

(0|a+afe|0) = Jv{q',,qy,E)v*{q^,E)e^piiiq'^-q,)k)dq'^dci^dE^ (116) 



The perturbation Hamiltonian, corresponding to the relaxation and decoherence processes have form 

Ahis (fc) = (fc) + (fc) , (117) 

where (see (f59|) ) 



(fc) ^ a/2 {I+S^ + I^S+) = (118) 



«4 



{u* {q^,E) ^+{qy,E^) + viq,^,E)C{qy,E+))expiiq,k) dEdqi_ + H.c. 



2(27r) 

The second term in Ea. (|117p 

A/ifj (fc) = -ahz {a+ak - (0|a+afc|0)) « -ahz (fc) ( S^S+ - ^) (119) 
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describes two-magnon interaction, which is similar to two-phonon interaction (Ref.|llj. sec. 3.4). This mecha- 
nism causes the modulation of nuclear spin resonance frequency without changing its state (adiabatic decoher- 

ence). It leads to the temperature depending decoherence rate, which is negligible small at (be — h^sBe ) ' 
when the system is very close to ground state. The zero-point electron spin oscillations are not involved in this 

(2) 

mechanism. Therefore, we will neglect next the contribution of terms A/ij^ (fc) and will use as the perturbation 
Hamiltonian the expression Ahis (fc) = A/ifj (fc). In this case, relaxation of transverse component of Bloch 
vector is accompanied by nuclear spin flopping (nonadiabatic decoherence). At the same time, the relaxation of 
longitudinal component of Bloch vector also occurs. Thus, these two processes may be considered here as one 
unified process of nuclear quantum state damping. 

We will now assume that interaction of nuclear spin, which is initially at coherent state (with nonzero 
no diagonal elements of density matrix pi{k,0)), with electron antiferromagnetic system in ground state, is 
switching on at the initial moment r = 0, when no perturbed density matrix is represented as direct product 
p{k,0)=pi {k,0)ps (0)=p7 (fc,0)|0)(0|. 

Let us go next to interaction representation for density matrix relatively to Hamiltonian ho (fc) = hs — 
{L0i{k)-a/2) h,: 

Pin {k, t) = exp {iho (k) t) p (fc, r) exp {-iho (fc) r) (120) 
and to the equation for density matrix of nucleus-electron system 

idpin {k, t) /dT = [Ahjs (fc, r) , Pin (fc, r)] , (121) 

where 

A his (k, t) = exp {iho (k) t) Ahjs (k) exp {—iho (k) t) = 

= a/2 exp {-i {loi (fc) - a/2) r) 1+ (r) + H.c. = (122) 

= ^ exp {-i {loi (fc) - a/2) r) /+ j [u* (q±, S) exp {-iE^t) ^+ {qy.E^) + 

+v (qj_, _E) exp {iE+r ) ^ {qy, E+)] exp (iqxk) dEdc[± + H.c. 
It is follows from Eq. (|12ip in the second order of perturbation theory 

idpin (fe, r) /dT « [Ahis (fc, r) , p (fc, 0)] - (123) 

r 

-I I [Ahis {T,k), [Ahis {k,T'),p{kM]dT'. 



Let us write the derivative of expression (|113|) with respect to time, by using Eq. (|123p and perform then the 
cyclic permutation under tracing. Finally, accounting the relation tr ([/^T, A/i/5 (fc, r)] ps (0)) = 0, we will find 

d [P- {k, t) exp {i {iui (fc) - a (1/2)) r)] /dr = 2 tr{rdpin {k, t) /dr} w 

r 

« -2 tr{/,7 / [Ahis (fc, r) , [Ahis {kr') , p (fc, 0)]] jdr' = (124) 



r 

^-2tr{J [[l,,Ahis{k,r)], Ahis {k,T')] p{k,0)}dT' . 



By defining the Bloch vector transverse component in the form 

P- (fc, t) = P- {k, 0) exp [~i {uji (fc) - a/2) r - (fc, t)] , (125) 

we will obtain: 

d [P- {k, t) exp (i {lui {k) - a/2) r)] /dr = (126) 

= -dr± (fc, t) /dr ■ P- (fc, t) exp {i {uji (fc) - a/2) r) , 

where ReFj_ (r) is decoherence decrement and ImF_L (r) is phase shift. 

We will represent the nuclear density matrix in the right part of Eq. ()124|) by expression 

PI (fc, 0) « 1/2{1 + 2P, (fc, 0) + P- (fc, 0) /+ + P+ (fc, 0) /- }|0) (0| |. (127) 
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Using in Eq. (1124|) the perturbation Hamiltonian ()118p . ignoring the insignificant factors 
exp (±i (luj (k) — a/2) r) and taking in the context of second order of perturbation theory in the left-hand side 
of Eq. ()126ll {k, t) P~ {k, 0), for the decoherence rate we will obtain 

T 

dRer±{k,T)/dT = RctTi J {0\[[l^,Ahjsik,T)], Ahjs {k,T')] I+\0)dT' = 



T 

Retr, / (0| [[47, {I^S- (r) + I' S+ (r))] , {l+S- (r') + I' S+ (r'))] /+|0)dr' = (128) 



4 

b 

} 

= -2Re j {0\S^ (r) 5+ (r') + S+ (r') (r) |0)dr'. 


where tr/ is partial trace over nuclear spin states. 

Let us restrict next again to low magnon excitation mode with energy E^. Further, let us insert week 
magnon damping {E^ E^ + is, E^ 3> s > 0) and make a set of rearrangements, which are similarly to 
that in Eqs.(A2.3),(A2.6). We will obtain 

Redr_L (fc, r) /dr « (129) 
^ ^ Re2 u* {c[±,E) u {q'^,qy,E) exp (— (ii?_ — s) (t' — r)) exp {i [q^ — q^) k) dq'^dq±dEdT' . 







(4^) 

Upon integrating over r', we will have 

Redr^ {k, t) /dr w (130) 



2 Re2 / u* (q^, S) u {q'^,qy, E) exp (z (g. - g^) k) ^ exp ((»g ^) ^) rfj^rfq^rf^;. 



(47r) 

By using the expression E (q^ ) w ■\/6p + q'^/12, we write 

dq^ = 27rq^dqi = 247ri; (q^) dE (q^) , < q_L < tt (131) 

and take notations A6fc = he — h — gk, £, ~ E — be- Upon integrating over E and q'^ (similarly to Eqs. 
([67]) - ([TO]) ), we will transform the expression (|130p for decoherence rate to the following form 



Redr^{k,T)/dT^^R^{Abk,T) = ^ I (^^Jl + bl + bc+^^ Y{^,+Abk,T) d^, (132) 



where 



and 



V|; + A6.,., = R.e '-"^'i-« + "''7'-l = ,133) 

- (n? + Abk) - s) 

(^ + Abk) sin ((^ + Abk) t) exp (-sr) ^ ^ 1 - cos {{£, + A6fc) r) exp (st) 
iC + Abkf+s^ ' (e + A6fc)Vs2 

1 — exp (— t) 



Ri_ (A6fc, r) = Rez ^ (^^Jl + b^ + be ~ Abk + i {t/r - s) 



dt, (134) 



with 



h = i-iAbk + s) T, t2^y~i\^^Jbl + n^/12-bc + Abkj+sjT. (135) 
Upon integrating over we have 



i?_L (A6fc,T) = Re{i Jl + + foe- - A6fc - is 



—I 
In- 



-iA6fe + s 
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Eiii-iAbk + s)T)^ El 



bl + n^/12~ be + Abk ]+s]t 



(136) 



where 



-l/Texp{-sT + iAbkT) |^exp \^\^i^^bl + tt'^/12 ~ &c ) t ) - 1 )}, 

OC 

El (x) = J ^^P(~^) rft, I arga;| < tt 



(137) 



is exponential integral. 

For the case of Ab^ 3> s > 0, omitting in dominator, we will have the explicit approximate expression 



R± (A6fe,r) 



l + bl + bc- A6* 



i {AbkT) - si Jbl + nyu -bc + Abk]T 



-1/r 



cos (AbkT) - cos I I Jfe^ + 7rVl2 - &c + ] t 



Abk{^ 



1 + 52 + 6c 



cos {AbkT) 
Ah 



+T ( si {AbkT) - si ( ( + 7r2/12 - 6, 




exp (— st) + 



^bl + nyi2 - be 



+ 7:^/12- be + Abk 
7r2/12 -bc + Abk) r 



(138) 



y%+^Tyu - be + Abk 



+ 



•c ■ 



Ab, 



■k T 



cxp (— st) > + 



+s < In 



ci {AbkT) — ci 



Tr2/12-&c + A&fc r 



exp (— st) > > 0, 



where si {x) = — J ^^dt, ci {x) = — / ^^dt are sine-integral and cosine-integral, and it is acounted, that 

X X 

lim [ci (ar) — ci (6t)] = In |. 

We notice (Fig. [7]) that oscillating part of the rate of decoherence tends fast (microseconds) to constant. 
Taking then r— >c» (r^l/s~ 10^), we will obtain for inverse time of decoherence (it corresponds in 
NMR science to transverce relaxation rate I/T2) 



3a2 

I/Td = iOE^cdV {k, 00) /dr = cj^— — i?^ (6fe, 00) 

ZTT 



(139) 



a2 ( + b''e + be- Abk) y/^^ + 7r2/12 - be ^^ ybl + nyU-be + Abk 
A6fcV&2^ + 7r2/12-6c + A6fc ^ 



A6fc 



For values a2 ~ 10^^, s ^ Afe^, loe/2tt ~ 10^^ Hz, decoherence time in sedconds is Tu ^ Abk/ s. 

It follows that decoherence time caused by one-magnon nonadiabatic processes near turning points is fast 
decreasing (for Abk < 10~^, s ~ 10~^ decoherence time may be less than milliseconds). Therefore, quantum 
operation, such as two-qubit operation, near these points should be carried out in the more short times. 

Between other mechanisms of decoherence it should be pointed out the adiabatic mechanism that is de- 
termined by magnetic interaction of nuclear spins-qubits with electron and nuclear spins of impurity atoms 
(Ref. [llj. sec. 5.4). The suppressing of this mechanism calls for appropriate cleaning of substrate from impurity 
atoms and using very low spin temperatures for nuclear spins. 

Note that character of decoherence rate essentially depends on the anisotropy of antiferromagnet (through 
value be) and from inhomogeneity of external field (trough parameter Abk)- 

The expression for frequency shift agrees at with Eg. ([55)1 for U {k, k) {s <C Abk): 



ImdF^ (fc,oo) / dT 



' 27r ^ 



+ be~ Abk In 



^bl + TTyi2~be + Abk 
Abk 



^2/12 -6c} 



(140) 
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Figure 7: A The r- dependence of Rj_ {Abk,T) (r = tw^;), for, 6^ = 1/4, s = 10"^ and Abk ^ be - bk = 
—3 • lO^'^, 1 • 10^"^, 3 • 10^"'. The case B corresponds to the increased values of R± (A6fe,r) only for Abk = 
1 • 10^3^ 3 • 10^3 in the range of 10"-* < r < 10"^. 
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Let us consider next the relaxation of Bloch vector's longitudinal component. Wc will write: 

BP, {k, t) Idr = 2 tr h^Opm {k, r) /Or = 

f T 

2 tr I 4, / [Ahis {k, t) , [Ahis {kr') , pi„ (fc, r')]]dr' I « (141) 



«-2Retr/ [[h,, Ahis {k,T)] , Ahjs {k,T')] pi {k,0)\0){0\dT' 

Taking {k,T) = P^ {k, 0) exp (— | (fc, r)) , we will obtain in the second order of the perturbation theory: 

dP, (fc, r) /dT « -dr,! (fc, r) /dr • (fc, 0) , (142) 
where the rate of the longitudinal relaxation is 

dr,! (/c,t) /dr = (dr,! (fc,T) /dr)* = 



a2 



2Retr/ y (0| [[I,, Ahis ik,T)] , Ahis{k,T')] /,|0)dr' = (143) 



2Retr, / (0| [[/., {I+S^ (r) + J- S+ (r))] , {l+S^ (r') + (r'))] 4|0)dr' 







4 

Similar to the foregoing calculation fEa. p28p ) we will have 



dr,! (fc,T) /dr ^ ^2Rc J {0\S+ (r) S~ (r') + 5" (r) 5+ (r') |0)dr' = (144) 



= Redr_L {k, t) /dr, 

whence it follows that for the considered mechanisms the rate of relaxation of longitudinal component is equal 
practically to the rate of relaxation of transverse component (decoherence rate). 

9 Decoherence of pair entanglement of quantum states in NMR 
quantum register 

The arbitrary state of pair removed spins-qubits in quantum register with zero Bloch vector values is described 
by the following reduced density matrix of nuclear spin system (a, P — x,y, z) (Ref . pjj , Sec. 2.5-2.7): 

p/(;,fc,T)= tr p/(l,...,m,...iV,r) = l/4{l + V4G„,;3(;,fc,T)(/j,„®/fe0)}. (145) 

m^Lk ^ — ^ 

Let us consider here the non-steady reduced density matrix of two nuclear state 

PI {I, fc, r) = 1/4 {1 + 4G,,, {I, fc, r) {Ii,Jk,z) + G+- {I, k, r) (/,-/+) + 0'^+ {I, k, r) (l+I^)}, (146) 
where diagonal and non-diagonal elements of the matrix are determined by magnitudes 

G+-- (Z, fc,r) = (G-'+(;,fc,T))*- 4 tr, (/+/-) pi{l,k,T). (147) 

{l,k,T) =4tr/ {Ii,Jk,z) PI (l,k,T), 

The Bloch vectors for this state Eg. piS)) are P {k, r) = 2 tr/ IkPi {I, fc, t) = P {I, t) = 0. 
In the interaction representation relative Hamiltonian 

ho {I, k)^hs^ {{uJi [l) - a/2) h (0 + (^/ (fc) - a/2) h {k)) (148) 
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the reduced density matrix has the form 



Pin (/, k,) = exp {iho (/, k) r) p (/, k, t) exp {-iho {I, k) r) . 



(149) 



Let there be the pure triplet entangled state of two removed spins I and k with zero total z-projection / = 1, 
M = 0, which belongs to the same sublattice, realized by the certain external action in the initial moment 
T = 0. It will be described by reducer state vector |1, 0) = ^1/2 ( | ti) + I it)) and density matrix: 



PI (/,fc,0) = |1,0)(1,0| 





110 

110 





1/4{1 - 4 + 2 (/,-/+) + 2 (I+Ife-)}, 



(150) 



with the following tensor components Gz,z {I, k, 0) = — 1, G+'^ (/, fc, 0) = G^'+ {I, k, 0) = 2. The concurrence of 
this entangled state has maximum value C = 1. 

Let us suppose now (as previously in Sec. 8), that the interaction of nuclear spins in ground coherent state 
pSOp with magnons is switching on at initial moment r = 0, when non-disturbed density matrix has the form 
of direct product: p{0) = pi {I, k, 0) ps (0) = pi {l,k, 0) \0){0\. 

We will write the equation for density matrix in the interaction representation pi^ {l,k,T): 



idpinil,k,T) /dr = [A/i/g (Z, fc, r) , pi„(?,fc,r)]. 



(151) 



The evolution of such two-qubit state in considered model is described by the following perturbation Hamil- 
tonian in interaction representation 



Ahis {I, fc, r) = [Ahis {I, t) + Ahis (fc, r)) = 

= exp {iho {I, k) r) {Ahis (0 + ^^/s {k)) exp {-ih^ {I, fc) t) 
where values Ahis {I, t) , Ahis {k, t) are determined by expression (|122p . 



(152) 



From Eq. p5ip for density matrix in the second theory of perturbation theory it is follows the equation 

idpin {I, k, t) Idr = \Ahis {I. fc, r) , pi„ {I, fc, 0) |0) (0|] - 



i / [A/i,s(Z,fc,r),[A/i,s(;,fc,T'),P/(Z,fc,0)|0)(0|]]dr'. 



(153) 



By using Eq. (jl53| ) and making the cyclic permutation, we will obtain the equations for the elements of 
reduced density matrix (|146p : 



dG,^, {I, fc, r) /dT = 4 tr/ {lLJk,z) 9pin {I, k, r) /Dj 



(154) 



-4tr/ 



and 



[[{Ii,zlk.z) , Ahis {I, fc, r)] , Ahis {I, k, r')] Pi {I, fc, 0) |0)dT', 



dG+^- {I, fc, r) /dr = 4tr (/j+^T) 9pi„ {I, fc, r) /di 



«-4tr, J (0| [[{l+I^),Ahis{l,k,T)] ,Ahis{l,k,T')] pi {l,k,0) \0) dr' . 



Taking into account Eg. pSOp . let us define the tensor longitudinal component in the form 

G,,, (Z,fc,r) = -exp[-r||(;,fc,T)] 
and the tensor transverse component in the form 

G+ - {I, fc, r) = 2 exp [i {lui {I) - lui (fc)) t~Tj_ {I, fc, r)] . 



(155) 



(156) 



(157) 



To obtain the expressions for the rates of pair relaxation of longitudinal and transverse components in the 
context of second order of perturbation theory, we will write: 



dGz,z {I, k, t) /dr = dPii {I, fc, t) /dr = dVJ^ {I, fc, t) /dr 



(158) 
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= -4 ti7 J (0| [[ilLJk,z) , ^his (l, k, r)] , Ahis (l, k, t')] (- (Ii^Jt^,) + 1/2 (/+/-) + 1/2 |0)dT' = 



= dFii (/, r) /dr + dT^ {k, r) /dr + dT^ {k, l-k,T) /dr 

and 

RedG+ - {I, k, t) /dr = -2Redrj^ {I, k, r) /dr = (159) 

r 

= -4Retr, J (0| [[{l^I^) , A/i,5 fc, r)] , fc, r')] (- (1^4..) + 1/2 (/+/,7) + 1/2 (/,-/+)) |0) dr' = 



= -2Redr_L {I, t) /dr - 2Redr_L {k, r) /dr - 2Redf _l {k, l-k.r) /dr, 
where the correlation part of longitudinal two spin relaxation rate is 

dFii (k, I — k,T) /dT = 

T 

= -4trz j {Q\[[Ii,,,Ahis{l.T)]h,.,Ahis{k,T')]{{l+I-) + {l-I+)) |0)dT'= (160) 



r 

= -^4Re j (0|[5+(r'),5rW] 10)^^' 



and the correlation part of transverse two spin decoherence rate is 

Redf _L [k, l~k,T) /dr = 

T 

= -2Retr/ ^ (0| [[/+, A/i/s (/, r)] A/^/s (fc, r')] (1;,./^) + (161) 



+ [[Ik,^his (fc,r)] I+,Ahis {l,r')] {iL^h,..) \0)dT' = 

T r 

= -^2Re / (0| [S- (r') , 5+ (r)] + [S^ (r') , S- (r)] |0) dr' = -^4Re / (0| [5+ (r') , (r)] |0) dr' . 



So, the correlation part of decoherence rate is equal to the correlation part of longitudinal relaxation rate. 

Let us rework next the expression (|16ip similarly as it was made in Sec. 6, accounting only low magnon 
excitation mode with energy : 

Redf_L (fc, l~k,T) /dr « 

« " ^ 2Re u* (q±, S) u {q'^,qy,E) exp (- (iE^ - s) (r' - r)) exp [i [q^k - q'J)] dEdq^dq'^dr' = 
(47r) J J 

o 

= -^2Re / u* (qx, i?) u (q^, g„ E) ^ - - s) r) ^. ^^^^ _ ^, ^^^q^^^, ^ (^g2) 

47r) J ^ [i-ty- — s) 



( V'rTfe^ + (q±) j - (-B(q^) - 6(i+fc)/2) sin(qj; (/ - fc)) + scos(9:r (^ - fc)) 



(47r)2y i?(q±) [ (i?(q±)- W)/2)' + -^' 

{E (q_L) - &(i+fc)/2) sin {{E (q^) - ?^(i+fc)/2) t - q^ {I - k)) 



{E{<i±)-\i+k)/2) +s 
3S ((g (q^) - 6(;+fc)/2) T-q^{l- k)) 
(^(q±)-60+fc)/2)Vs^ 



exp (— sr) > dqj^. 



27r 

Taking next in to account that J sin [qj^ cos ip (I — k)] dip = and introducing again the variable £, — E (q_L) — 



be, we will write E (qj.) — b(i+k)/2 — E (q±) — hk ~ g {I ~ k) /2 ~ ^ + Abk — g{l ^ k) /2 and for the correration 
part of decoherence rate we obtain (Fig. [5]) 

Redr_L (fc, l-k,T) /dr « —Ri_ {Abk, l-k,T)= (163) 

2tt 
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Figure 8: The r— dependence for the correlation part of dimensionless decoherence rate R± {Abk,l — k,T) for 
values A6fe = 3.10"^, s = 10"^, / - fc = 200 (away frow turning point), l-k = 299, 300, 301 (in the vicinity of 
the turning point). 



/ 



(^^Jl + bl + bc+^^ Y{^ + Abk-g{l-k)/2,T) Jo Uu [{be + 0^ ~ b^] {I - k)\ 





The decoherence rates of entangled qubit pair are due to decoherence of one spin states I, k and also 
to the correlation between nuclear spins /, k. The diagonal and non-diagonal elements of density matrix 
1 + Gz,z {I, k, t) and G^'" (Z, /c, r) are decreased with full rates RedFn [l, k, t) jdr and RedFx (^ fc, t) /dr. Note, 
that the asymptotic value for correlation part of decoherence rate (Fig.[S]) may change sign, if Ab^—g (I — fc) /2 < 
0, that is for states after "turning point" (Eg. ((75)) ). 

In the case being considered the concurrence for entangled two-qubit state can be obtained by using the 
Wootters formula (Ref.p^]) 

C (/, fc, r) = 1/2 max{|G+'- (?, fc, r) | - (1 + G,,, (/, fc, r)) ; 0}. (164) 

Taking into account that ReFj^(fc,T) — F||(fc,r), |G+'^ (Z, fc, r) | « 2 exp (— RcF_l (/, fc, t)), Gz,z {l,k,T) fa 
— exp (-F|| (/,fc,T)) and also Eqs p^ . (fTST]) . (fTsg]) . we will obtain 

C{l,k,T) = l/2{3exp(^-ReFi(?,r)-ReF_L(fc,T)-Refx(fc,Z-fc,r)) -1}. (165) 
For the logarithmic damping rate of concurrence we will then write 

d In G (/, fc, t) /dr = -'d/2d f RcF^ (^ r) + ReFj^ (fc, r) + Ref x (fc, I - fc, t)) dr. = (166) 



= -^iRi- (Afo;, r) + (A6,, r) + {Abk,l - fc, r)) . 

ZTT 

For A6fe > parameter A6/ = A&/j — g (I — k) changes the sign and parameter in denominator may not be 
omitted. Using next Eq. (|133p we will write 



R^{Abi,T)= J (^^l + bl + bc+^^ Y{i + Abk~g{l~k),T)d^ (167) 
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5.6 

/^(A6,,/-/c,Tj 




Figure 9: The r— dependence of (A5fc, I — k,T) for values Abk = 3.10 ^, s = 10 ^, I — k = 200 (away frow 
turning point), I — k — 299, 300, 301 (in the vicinity of the turning point). 



FinaUy, for the value of concurrence-damping rate in the context of second order of permutation theory we 
will obtain 



where (Fig. O 



dC {I, k, t) /dT = -—Rf {Abk, l-k,T), 
Ztt 



Rf {Abk,l -k,T)^ R^ {Abk,T) + i?i (A6/, r) + i?^ {Abk, I - k,T) = 



1 + bl + bc + n ■ [Y {^ + Abk,T) + Y {^ + Abk - g {I - k) ,t) 



+ Y{i + Abk-g{l-k)/2, t) Jo (^^12 [{be + 0^ - {I ~ k) 

The velocity damping of concurrence tends to constant for r — > oo: 

Rf {Abk,l-k,(x) = 

' 6c + e 



dt 



{^ + Abkf {^ + Abk-g{l^k)f 



Jo Jl2 



{be + if -hi {l-k) 



{i + Abk-g{l-k)/2Y 



di > 0. 



(168) 



(169) 



(170) 
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10 Appendixes 

Al. Determination of the coefficient A in Eq. (l5! 

Substituting Eqs. (|50p . (|5ip in condition p2p . we will obtain the equation for determination of the quantity A: 

J cxpS^-'-{{E-E')q,)yM{E',q^)dq., + 0{g) = S{E-E'), (Al.l) 



2EA^ 



^nn^+E 

where 



M(i;,qx) = 5I^^P| "J I / E{q^)dq,-EQ, 



771— — oo 



^ exp<j--| / Eicii_) dq,-E2nm\\. (A1.2) 



To achieve the convergence of summing over m in Eq.(Al/2), we will introduce a weak damping: 

iEm ^ iEm - i^\m\, ^ +0. (Al.3) 

Using Poisson summing formula (Ref.|23j. Eqs. 4.8-9) and performing the rearrangement of this sum to 
integral over continuous variable 27rm — s- gfi, we obtain the following asymptotical expression 

oo oo / \ 

M {E,c[±_) = Um^ g J exp{-n J E (q^) dq^ - E fi - p^i\ - lylfi]} dfi 



^ lim f exp{-j {E (q^) - E) pi/g ~ ,y\pi\/g} dfi + O (.g)} = g[S{E-E (q^)) + O (g)] . (Al.4) 
27r u^+a J_r^ 



The last result ignores here the fast oscillating terms with p ^ 
On the other hand, we may write Eq.(Al.4) in the form 

M 



'I(E,qA= lim V exp{-- ( / E (q±) dq^ - E2TTm\ - -v2TT\m\} + O (g) = 
mt^oo 9\ J 9 

oo { ■ I 9^+2'^™ \ 

= 1 + 2Re lim V exp < -- / E iqA dq^ ~ E2TTm - v2'Km > + O (.g) w 1 + O (g) , (Al.5) 
--^t. \ ^ \ I )\ ■ 

where the oscillating sum is omitted. 

Comparison of Eq.(Al.4) with Eq.(A1.5) allows making the conclusion: 

lim g5 (E - E (q^)) ^ I (A1.6) 

To achieve the convergence of integrating over qx in Eq.(Al.l) we introduce again a weak damping: 

lE'qx^lE'qx-s\qxl s +0. (Al.7) 

Let us next extend the integration limits over "extended" variable rj — qx/g to±oo. Then the expression in 
the left part of Eq.(Al.l) takes the form 



1"°° i 

lim / exp{'-iE~E')qx~s\qx\}-MiE',q^)dqx 



= lim g exp{i{E-E')r]-s\Tj\}-{'^ + 0{g))dT] = 2ngS{E''E'){l + 0{g)), (Al.8) 

"^+0 J-oo 

where result (Al.5) and the relation 

/oo 
exp [i {E ~ E')ri~ s\r]W di] = 2TrS {E - E') (Al.9) 
-oo 

was used. 

For determination of value of A^ we will now obtain from Eqs. (Al.l), (Al.8), the equation: 

AnqEA^ 
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A2. The calculation of indirect nuclear spin interaction 

After commutator transformation in Eq. ()65p we will find 

hii {k,l) « 



« -mV4{/+ (/c) /- (0 J (0| [S'-(r,fc), S'+(T + r',0] |0) exp(sr')rfT'+ 

— oo 



+I+{l)I-{k) ( (0| [S'-(r,0, 5'+(r + r',fc)] |0) exp(sT')dT'} +H.c.- (A2.1) 



u 

-zaV4 5^ / {l/2{Q\[S-{T,j), S+{t + t',j)\\Q)- 

j=k,l_ 



(2t 



+h (j) {t,j) S+ (r + t',j) + S+{t + t\j) (t, j) |0) exp (sr')}dT' + H.c, 

where it was accounted that 

U) U) - 1/2 ± (j) , [I^ U) , (j)] = ±2/. (j) , (A2.2) 

(j) ij) - -I^ ij) h (j) - ±l/2/± (j) , [/. 0') , /± (j)] = ±/± (j) . 
The mean value of commutator of electron spin operators takes the form 

(0| [S'{T,k), S+{T + r',l)] |0) = 

^ [ {u*{q^,E)u{c^^,E') (0| [t{qy,E^),^{q'y,E'_)] \0) exp {-iE'_t' + sr') + 
in) J 

+v {q^,E) V* (ql, i?') (0| (qy, E+) , ^ {q'y,E'_^)] |0) exp {iE'_^T' + st')}- (A2.3) 
• exp [i (qxk — q'J)] dEdE' dq±dq'j^. 
We will use next the mean value on ground state of commutator of magnon operators (I2ip . written for r = 0: 

(0| (qy, E) , e {qy,E')] I 0) = (0|e (qy, E) ^ {q'y.E') | 0) = 5 - q'y) S {E - E') (A2.4) 

and thus we will have 

(0| [5-(r,fc), 5+(r + r',0] |0) = 

{tf I '''''' (-*^-^' + *^')+ (A2.5) 

+v {q±,E) V* {q'^, qy, E) exp {iE+t + st')} ■ exp [i ((7a;/c - q'J)] dEdq±dq'^. 
We will find then for mean value of the electron spin operators product the expression: 

(0|5- (r, j) S+ (r + t',j) + 5+ (r + t', j) {t,j) |0) = 

= 7A2 / (q^' ^) ^ il'^^ly^ ^) (~iS-r' + st') + (A2.6) 
(27r) J 

+v{q^,E) v*{q'^,qy,E)exp{iE+T' + ST')}exp{i{q^-q'Jj) dEdq^dq'^ 
and Eq.(A2.1) takes the form: 

hii (kj) = 

^ 

= -*T7T^/ / {I^ik) I- il)[-u*{q±,E)uiq',,qy,E) exp i-tE.T' + ST') + 
4 (27rJ J J 

—00 

+v* iq±,E) V {q'^,qy, E) exp {iE+t' + st')] exp [i {q^k - g^^] + 

+/+ (/) /- (fc) (q±, £;) M £;) exp i-iE^T' + st') + (A2.7) 

+v* {q±,E) V {q'^,qy,E) exp {iE+r' + sr')] exp [i [qj - g^fcjJjrfE'dqj^dg^dr'- 
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4(2^)' 



J2 {1/2 [-U* (qx, u {q',,qy,E) exp {-iE^t' + st') 



+v {q±,E) V* {q'^, qy, E) exp {iEj^t' + st')] + 

+4 [u* (q_L, -E) u {q'^,qy,E) exp {-iE.r' + st') + 

+w (q±, w* {q'^,qy,E) exp (iS+r' + sr')] exp [i (q^ - q'^) j]}dEdq_i_dq'^dT' + H.c. 

After integrating of Eq.(A2.7) over r', we will obtain the expression for the effective two nuclear spins 
Hamiltonian for atoms k and I among the same sublattice A as a correction of the second order perturbation 
theory to the energy of antiferromagnet ground state: 



hii {k,l) = 

= - E {[^^ (j) - a/2 - 14- (j)] 4 + [/ ij,j) 12} - U {k, I) [I- (fc) /+ (0 + /+ (fc) /- (0] 
j=kj 



(A2. 



where the expressions for indirect interaction between two separated nuclear spins, belonging to common sub- 
lattice and generalizing the known Nakamura's expression for antiferromagnet in homogeneous fields [7j and for 
correction for nuclear spin resonance frequency have the forms 



2 (27r)' 



rRc 



U{k,l) = 

u* {qa:,qy,E)u{q'^,qy,E) v {qa;,qy, E) v* {q'^,qy,E) 



E^ 



IS 



E. 



IS 



and 



2(27r)' 



rRe 



w U) = 

u* {qx,qy,E)u{q'^,qy,E) v {q^,qy,E)v* (q'^,qy,E) 



+ is 



E. 



IS 



(A2.9) 

exp [i {q,j:k ~ q'J)] dEdq±dq'^, 

(A2.10) 

exp [i {q^ - q'^) j] dEdqj_dq'^. 



11 Conclusion 

In this paper it was considered, as a quantum register, one-dimensional chain of magnetic atoms with nuclear 
spins, which is placed by regular way in thin plate of easy-axis 3D antiferromagnet. 

When the external magnetic field is directed along the easy axis normally to the plane of the plate and has a 
constant gradient along the nuclear spin chain, antiferromagnet spin Hamiltonian in the case of inhomogeneous 
external field in spin- wave approximation was obtained. If the field has a weak gradient, asymptotic expression 
for coefficients of unitary transformations to diagonal form of this Hamiltonian was found. 

The expression for indirect inter-spin coupling that is due to hypcrfinc nuclear-electron coupling in atoms 
and spin- wave propagation in antiferromagnet was evaluated and it was shown that owing to gradient of external 
field significantly change. So, in the case that the value of local field in the middle point position for two spins 
coincides which critical field for homogeneous phase transition, the indirect interaction fast grows and, even 
so, it takes the oscillating character of distance dependence. We have denoted the corresponding points as 
"turning points" . In conditions of homogeneous AFR the indirect interaction has additional turning points, 
which determined by frequency and microwave power. 

The nonadiabatic mechanism of nuclear spin statesdecoherence in quantum register caused by interaction of 
nuclear spins with magnon excitations was considered and the calculations of one qubit and two qubit relaxation 
and decoherence rates were made. 

It turns out that the character of decoherence processes essentially depends on antiferromagnet anisotropy 
(parameter be) and on inhomogeneity of external field (parameter g). As this takes place, the temperature 
whereby the thermal magnon excitations are excluded and the two-magnon spin-lattice relaxation is especially 
suppressed, should be defined by values T <C Tc {1 — b/bc) , ?b = h'fsBE/kBBc- As an example we have 
also considered decoherenc of pair qubits maximally entanglement state and have calculated the concurrence 
damping rate. 

Due to the availability of inhomogeneous external field, it becomes possible to control not only the individ- 
ual nuclear spin resonance frequency but also the interaction between spatially separated spins in large-scale 
quantum register without resort to controlling many gate systems. The switching of interaction between far 
removed spin-qubit, required for two-qubit operation, in considered register may be performed by the tuning of 
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qubit state to turning point state, where indirect interqubit interaction has a large value. The external magnetic 
field, its gradient and microwave power may play the role of control parameters. 

For the examination of realization problems of considered model for large scaled quantum computer, as 
in the case of any other model, one should refer to five necessary fundamental requirements given in paper 
[24j . One of the most difficult problems in the realization of quantum computer is the problem of quantum 
state initialization. For its solution, it is believed to use the method of dynamical polarization of nuclear 
spin in antiferromagnet for homogeneous field condition in the manner similar to the method suggested for 
semiconductors in book '11] and the scheme of optical pumping, as the one discussed in paper |25| . The second 
difficult problem is the realization of readout processes, for this purpose an ensemble approach scheme like the 
one considered in book [TT] may be suggested. 
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